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Abstract. Numerical integration of the Landau-Lifshitz-Gilbert equation 
with thermal fluctuations is used to study the dynamic response of single- 
domain nanomagnets to rapid changes in the applied magnetic field. The 
simulation can resolve magnetization patterns within nanomagnets and uses 
the Fast Multipole method to calculate dipole-dipole interactions efficiently 
The thermal fluctuations play an essential part in the reversal process when- 
ever the applied field is less than the zero-temperature coercive field. In this 
situation pillar-shaped nanomagnets are found to reverse through a local 
curling mode that involves the formation and propagation of a domain wall. 
Tapering the ends of the pillars to reduce pole-avoidance effects changes the 
energies involved but not the fundamental process. The statistical distri- 
bution of switching times is well described by the independent nucleation 
and subsequent growth of regions of reversed magnetization at both ends 
of the pillar. 



Magnetic nanoparticles are important components of nanotechnology. 
Previously Ising models, which describe only highly anisotropic materi- 
als, have been used to understand thermally activated switching in single- 
domain magnets [I]. Single-domain magnets with dimensions of only a few 
nanometers are being manufactured and measured, for instance, using scan- 
ning microscopy techniques [2]. We have chosen these nanomagnets as a 
specific system in which to investigate magnetization switching with more 
realistic computational models of magnetic materials. 

Here magnetic materials are modeled by position-dependent magneti- 
zation density vectors M(r) with fixed length M s , that precess around the 
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Figure 1. Hysteresis loops of period 4 ns for fiat-ended nanoscale pillars described in 
the text. One is at T = OK (solid) and two at T = 100 K (dotted). 



local field H(r) according to the Landau-Lifshitz-Gilbert equation [3, 4] 

H(r) - ^-M(r) x H(r)] , (1) 

where 70 = 1.76 x 10 7 Hz/Oe is the electron gyromagnetic ratio and a 
is a phenomenological damping parameter. The value a = 0.1 was cho- 
sen to give an underdamped system. The saturation magnetization M s = 
1700 emu/cm 3 and the exchange length £ x = 2.6 nm were chosen to match 
those of bulk iron. Numerical integration of Eq. (1) was carried out using a 
finite-differencing scheme with Ar = 1.5 nm and At = 5 x 10 -5 ns. The local 
field includes exchange interactions, dipole-dipole interactions, and thermal 
fluctuations. Details of the numerical approach will appear in Ref. [5]. 

Here we consider nanomagnets with square cross sections of side 9 nm. 
One geometry has flat ends, and is 150 nm along the z-direction [6]. A 
simulated hysteresis loop at T = K with a period of 4 ns for this system is 
shown in Fig. 1. There is no crystalline anisotropy included in the model, 
but the shape anisotropy is quite strong, giving a coercive field of about 
1875 Oe [5]. Even for these extremely fast loops, for a significant amount 
of time the magnet experiences a field just below the coercive field, and 
thermal fluctuations can carry the magnetization through the configuration 
associated with the energy barrier. Some hysteresis loops simulated at T = 
100 K are also shown in Fig. 1 to illustrate the extent of this effect. 

A more dramatic way to see the effects of thermal fluctuations on the 
reversal of the magnetization is to hold the magnitude of the applied field 
fixed just below the coercive field. This traps the magnetization in a shallow 
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Figure 2. Snapshots of the ^-component of the magnetization, M z , at equally spaced 
times for (a): pillars with square ends and — 1800 Oe (b): pillars tapered at both 

ends and ff ap p = 1975 Oe. All are identical cutaway views with only half of the pillar 
shown. Light shades indicate magnetization pointing in the unfavorable direction while 
the darkest shades indicate magnetization pointing in the favorable direction. Both sim- 
ulations are for T = 20 K. 



metastable energy well, and thermal fluctuations become an essential part 
of the reversal process. Snapshots of the z-component of the magnetization 
separated by 0.05ns are shown in Fig. 2(a). The end caps associated with 
pole avoidance become strong in an antiparallel field, and thermal fluctu- 
ations affect their volume. The reversal starts when thermal fluctuations 
make one of these regions supercritical [6], after which it propagates at a 
constant rate towards the other end of the magnet. Tapering the ends of the 
magnet reduces the pole-avoidance effects, and hence affects the end caps. 
However, the switching is still characterized by nucleation at the ends; an 
example is shown in Fig. 2(b) for pillars 227 nm tall. The major effect we 
have observed of the tapering is an increase in the coercive field. 

A simple model based on independent nucleation at the ends, followed 
by constant growth, allows calculation of the probability of not switching 
before time t to be [5, 6] 
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Figure 3. Probability of not switching, P no t(t), for 100 simulations of the flat-ended 
pillars at H app — 1800 Oe and T — 20 K, along with a fit to Eq. (2). The inset shows 
similar results for 22 trials with the tapered-end pillars at -ff app = 1975 Oe and T = 20 K. 

where to is the earliest time at which a pillar can switch (when both ends 
start growing at i=0) and T is the nucleation rate for each end. This form 
is compared to the simulation results for both pillars with flat and tapered 
ends in Fig. 3. The biggest obstacle to testing Eq. (2) is generating sufficient 
statistics for the longest and shortest switching times [5, 6]. 

In summary, thermal fluctuations and dipole-dipole interactions have 
been included in a three-dimensional model of single-domain nanoscale 
magnetic pillars. Simulated dynamics using the Landau-Lifshitz-Gilbert 
equation with thermal noise show that magnetization switching in the pil- 
lars occurs as the end caps nucleate, grow, and eventually coalesce. A simple 
model based on the thermally-activated nucleation of end caps for subcoer- 
cive fields gives a reasonable description of the simulation results for P no t(t). 
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